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Abstract. We consider an array of lb Cooper Pair Boxes, each of which is coupled 
to a superconducting reservoir by a capacitive tunnel junction. We discuss two effects 
that probe not just the quantum nature of the islands, but also of the superconducting 
reservoir coupled to them. These are analogues to the well-known quantum optical 
effects 'superradiance,' and 'revival.' When revival is extended to multiple systems, 
we find that 'entanglement revival' can also be observed. In order to study the above 
effects, we utilise a highly simplified model for these systems in which all the single- 
electron energy eigenvalues are set to be the same (the strong coupling limit), as 
are the charging energies of the Cooper Pair Boxes, allowing the whole system to be 
represented by two coupled quantum spins, one finite, which represents the array of 
boxes, and one representing the reservoir, which we consider in the limit of infinite size. 
Although this simplification is drastic, the model retains the main features necessary 
to capture the phenomena of interest. Given the progress in superconducting box 
experiments over recent years, it is possible that experiments to investigate both of 
these interesting quantum coherent phenomena could be performed in the forseeable 
future. 
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1. Introduction 

Recently there has been much interest, botli theoretical Pl Ej and experimental 
mis El, in superconducting islands coupled to superconducting reservoirs through 
Josephson tunnel barriers. Such structures can be manufactured with energy levels 
such that the islands approximate two - level systems, and may have applications as 
qubits [7]. In particular, work has focussed on demonstrating the quantum mechanical 
nature of these systems relevant to quantum computing. A full demonstration of many- 
qubit quantum algorithms in these systems may be difficult at present, and it may be 
helpful to find more general indicators of the quantum coherent nature of these systems 
without the need to, for example, measure a Bell inequality. An ability to demonstrate 
coherence is a prerequisite for any quantum computer, and so could function as a helpful 
first step. In particular, revival demonstrates coherence on a timescale which is long 
compared to the simple oscillation period. This demonstrates that the decoherence time 
is long enough for multiple gate operations to be performed. 

In this paper, we discuss two effects that probe not just the quantum nature of 
the islands, but also of the superconducting reservoir coupled to them. The first is an 
analogue to superradiance in quantum optics 0, in which the intensity of the radiation 
emitted from a collection of atoms is proportional to the square of the number of atoms 
rather than linearly proportional. Arrays of Josephson junctions coupled through a 
cavity [HI have also demonstrated enhanced emission of radiation (Barbara et al jlUj). 
Here we consider the current emitted from an array of Cooper Pair Boxes coupled 
through a superconducting reservoir, which we call the super- Josephson effect. 

The second effect is the phenomenon known in quantum optics as revival. In this, 
the coherent superposition properties of an initial quantum state decay through coupling 
to the quantum degrees of freedom of the reservoir, and then 'revive' at a later time. 

When revival is extended to two or more two-level systems we discovered a new 
phenomenon: 'entanglement revival.' 

In order to study the above effects, we utilise a highly simplified model for the 
corresponding systems of superconductors. Although the principle simplification, which 
is to set all the single-electron energy eigenvalues to be the same, is drastic, the model 
retains the main physical features necessary to capture the phenomena of interest. For 
example, the essential feature needed to observe revival is the discreteness of the Cooper 
Pairs in the small superconductors and the reservoir. 

2. The Model 

2.1. Motivation 

The Josephson effect is an intrinsically quantum mechanical phenomenon, as is 
superconductivity itself. However, in the BCS approximation for a bulk superconductor 
the phase of the order parameter is usually treated as a classical variable ^T] , in the sense 
that fluctuations about it are considered negligible. When the size of the superconductor 
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is reduced, these fluctuations become relevant, and the mean-field treatment is no 
longer valid. Such deviations from mean field have been observed in experiments on 
superconducting nanoparticles by Ralph, Black and Tinkham (RBT) ^2]- Subsequent 
work in the field is well reviewed in |T^. 

In what follows we shall be concerned with the same fiuctuations about the mean 
field theory, albeit for slightly larger superconducting samples, ~ lOOnm, than in the 
RBT experiment [T^ . 

In general, to go beyond the BCS approximation is quite difficult even for a 
single bulk superconductor. Although an exact solution exists for a finite size sample 
[m ITK] , this returns a set of coupled equations which rapidly becomes intractable for a 
superconductor with more than a few levels. 

In our case the difficulty will be compounded by the fact that we wish to treat an 
array of small superconductors individually coupled to a large one. Thus, to render the 
problem tractable we shall introduce a highly simplified description of each component 
of our system. Nevertheless, we shall endeavour to retain sufficient realism to capture 
some generic features of a number of surprising novel phenomena. These arise when we 
study the collective quantum states of an array of small superconductors which are not 
coupled to each other but individually coupled to a common large superconductor; a 
Cooper Pair reservoir. The model will be introduced in subsections 12.21 and 12.31 The 
two new phenomena whose description is the principle aim of this paper, namely the 
superconducting analogues of the quantum optical phenomena of 'superradiance' and 
'revival,' will be treated in sections El and ^ respectively. 



2.2. The Cooper Pair Boxes 

The system we shall study is depicted, schematically, in Figure [T] It consists of an 
uncoupled array of small superconductors, which we shall refer to as Cooper Pair Boxes, 
each of which is coupled through tunnel junctions to the same bulk superconductor which 
we shall call the reservoir. The number of Cooper Pair Boxes will be denoted by 
The hamiltonian of the system is: 

H = Hb + Hr + HT (1) 

where the individual hamiltonians refer to the array of Cooper Pair boxes, the reservoir, 
and the tunneling between the two respectively. The first term is the hamiltonian of 
the boxes: 

H, = j:Efat (2) 

Each Cooper Pair Box is capacitively gated in such a way that it can be described 
as a two-level system where the two states are zero or one excess Cooper Pair on the 
box, with Ef^ the charging energy of each. In short, they are the objects investigated 
in the experiments of Nakamura et. al. |S1 E] • 
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Figure 1. An array of It superconducting islands, or Cooper Pair Boxes, that are 
individually coupled to a superconducting reservoir, r, through capacitive Josephson 
tunnel junctions (with no inter-box coupling). 



The second term is the well known BCS hamiltonian (see section . and 
the final term: 

Ht = -Y. {^tc-kiCk^ + fTr4^clfc^) (3) 

describes the tunnelling of Cooper Pairs from each of the boxes to the reservoir (this 
form can be derived from single-electron tunnelling ^Hl with Tj ^ proportional to the 
square of the single electron tunnelling energy, tfj?). The charging energy of each box can 
effectively be tuned by applying a gate voltage to that box, and a tuneable tunnelling 
energy can be achieved by adjusting the fiux through a system of two tunnel junctions 
in parallel If the charging energies Ef^ and tunnelling elements Ti^k of each box 
are all adjusted to be the same, then the charging energy for all the boxes acts like the 
z-component of a quantum spin with J = lb/2. Similarly the sum over all the raising 
(lowering) operators acts like a large-spin raising (lowering) operator. 



i=0 

t^t^S^ (4) 

i=0 

If we consider transitions between the states of the whole system of boxes caused 
by only these operators, then the symmetrised number states form a complete set. The 
eigenstates of represent states where a given number of the boxes are occupied, with 
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m = Nh — lh/2. We now have a description of the Cooper-Pair boxes in terms of a single 
quantum spin: 

Tj rpch qZ 

Ht=Y. -Tk{St c.kiCk^ + ^6~4|clfc^) (5) 

k 

2.3. The Cooper Pair Reservoir 

The quantum properties of the reservoir will also be of interest and hence we must 
simplify its description as well. To do this we start with the BCS Hamiltonian written 
in terms of the Nambu spins [T7], (erf = \{c\^Ck^ + cl;,|C_fcj, ~ 1); cr^ = c\^c^_^,y): 



k k=kp—uic 

(6) 

where is a generic label (not necessarily the wavevector) of the single-electron 
eigenstates and the pairing field, for an electron-electron coupling constant V, is given 
by Afc = V{a'^) = V {c^kiCki) ■ As usual, this has to be determined self-consistently. 

We shall now make the approximation that the pairing fields, Afc,A^, are 
independent of k, and that the term 2(6^ — fi) is also the same for all k. Remembering 
that the interaction only occurs in a small region around the Fermi energy, ep ± Uc, 
the approximation that all the are equal is equivalent to taking the limit Uc V 
Appendix A| ). Clearly this is a drastic, strong coupling, approximation P^IT^ I^. 



see 



but we will demonstrate that it nevertheless retains some important features of the 
superconducting state. The gain in making this approximation is that now the Nambu 
spins can be treated as a single large spin. Further discussion of the justification of 
this model when the superconductor is finite and its relation to the Richardson solution 
can be found in ^HI- Note however, that Yuzbashyan et. al. jTHI are concerned only 
with isolated islands with a constant number of Cooper Pairs and therefore discard the 
diagonal term {K^). 

This representation of a sum over Nambu spins as a single spin was used by Lee 
and Scully in 1971 but that paper discarded the commutator [S^^f^, S~f^] = 25"^^. 
By retaining it, we can deal with the discreteness of the Cooper Pairs. The size of the 
spin is given by 5*^ = /r /2 where Ir is the number of levels within the cutoff region, which 
is proportional to the volume and will be taken to infinity at the end of the calculation. 

In the above limit the hamiltonian ^ becomes: 

Hr = 2^FrS^ - S^A - S;A* 

(7) 

Where {epr — A*) is written as ^pr for convenience. 
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Before treating the coupled system of boxes plus reservoir, we investigate the 
behaviour of the reservoir alone, in order to demonstrate that our description still makes 
sense even after the above approximation. 

The hamiltonian of the reservoir ((Tj) is composed of the , S'^ and 5*7 operators, 
and so can be written as the component of spin along some other direction h. 

Thus we attempt to diagonalise Hr by defining S^f^ by •yrSrh = -^r- Then the 
commutation relations [S*^^, 5*;^^] = S^^f^ and [S^^f^, 5*"^] = 25*^- lead us to the raising 
and lowering operators in this direction: 



\d\ = \A\/Ef 
-dA 

6 = 

— 2,Ep 

^ -dA* 

(8) 

Where 7^ = Ep = ^JQ^^ + I^P; i-^- quasiparticle energy evaluated at the Fermi 
energy. 

The ground state of the hamiltonian is, obviously, the m = —lr/2 eigenstate in the 
new basis. In the old basis, this is the spin coherent state [22] : 



1 ((y*S+)^ 



which can be rewritten as: 



\c^)= I \ , na+«V)io) (10) 
\/i + i"i k 

As can be readily ascertained using S^^f\a) = —lr/2 \a), choosing a = ^•^''^'^^ 
ensures that this is the ground state. 

In the second expression for \a) we have explicitly rewritten the coherent state as 
a product over k. In this form it is clear that this state is identical to the BCS state in 
the limit <^ V, when Uk = Vk = IV2. The pairing parameter A and the chemical 
potential /i are determined self consistently. In the conventional BCS theory [TT] : 
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For our simplified spin model, the fact that the are equal allows us to evaluate 
the sums in flllj) . giving two coupled algebraic relations instead of integral equations. 
Thus we can evaluate A and fi exactly. We find: 



|A|2 = V^Nrilr - Nr) (12) 
eFr-f^ = V{lr/2 - Nr) (13) 

In order to check the validity of these relations, we solve the gap equations for a 
general BCS case, and then take the limit where all the levels have the same energy. To 
take this limit, we consider a top-hat density of states centred around ep with width w 
and height l/w. This allows the limit where the bandwidth w — > to be taken whilst 



keeping a constant number of levels, /. As shown in Appendix A at zero temperature 
the integrals can be calculated exactly. We find: 



A-(..)^(|±f-l).(f)^(|^-.) 

iV. = ^ (I - ^i^Fr + + + ^i^Fr " |)^ + A^) (15) 

With X = e^'-r and n = -r^. These two equations can be solved to give A and ^pr- 
When the limit V/w oo (the strong coupling limit) is taken, we find that we obtain 
the same results as given by the spin model, (fT^ IT^ . 

Evidently, the fact that this spin representation reproduces the BCS results, albeit 
only in the strong coupling limit, lends credibility to our simplified model and allows 
us to consider the situation we wish to describe, namely that of a series of Cooper Pair 
Boxes coupled individually to a single superconducting reservoir. The problem is then 
that of two quantum spins, one finite, which represents the array of Cooper Pair Boxes 
(labelled b) and one for the reservoir (labelled r) which will consider in the limit lr oo 
at the end of the calculation. In short, we shall consider the hamiltonian: 



H = E^S^ 

+ 2iFrS^ - S^A - S;A* 

- T(5+5- + S^S,) (16) 



3. The Super - Josephson Effect 

The first effect we consider is an analogy of superradiance (see e.g. In quantum 

optics, this refers to the fact that lb atoms, placed in a particular state and coupled to a 
mode of the electromagnetic field will emit radiation proportional to the square of the 
number of atoms. The initial state is one of the form: 



|7/>(0)) = {|000..111..) + I001..110..) + ■ ■ ■ + I111..000)} 



(17) 
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where and 1 indicate that an atom (identified by the place of these numbers in the array 
which specifies the state ?t,2, ?t,3 . . .)) is in its ground or excited state. Clearly, this 
state is an eigenstate of the operator for the total number of excited states (S^), but not 
of the individual af operators. Note that it is the symmetrised state with a total number 
of excited states Nb. Evidently, these are exactly the states described by the model in 
section 12.21 It is also interesting from the quantum information theoretic point of view 
that these are highly entangled states. These states could, in principle, be generated by 
tuneable inter-box coupling [2 . For a carefull discussion of superradiance the reader is 
referred to the very readable account of Eberly [2Sj- Note also that Josephson junction 
arrays placed electromagnetic cavities can emit radiation with similar dependence 

uni. 

In this paper, rather than electromagnetic radiation emitted from atoms (or 
Josephson junctions) coupled through a cavity, we consider the current that is emitted 
from an array of Cooper Pair Boxes coupled through a common superconducting 
reservoir. We consider three cases. The first is when the system is in an eigenstate of 
both the number operator of the large superconductor, and the operator measuring the 
total number of Cooper Pairs on the boxes. In this case, whilst our model incorporates 
the fact that the electrons are paired, there is no phase coherence between these pairs. 

This is an artificial case if we think of the boxes coupled to a bulk superconductor, 
but instead would represent many boxes coupled to an island with only a few levels. If 
both the box and the reservoir are finite then, of course, the order parameter | A| is zero. 
The square of the order parameter |Ap, however, is non-zero, indicating the presence 
of pairing between the electrons. It is, however, the closest to superradiance, where the 
atoms are in a symmetrised number eigenstate, and the field contains a given number 
of photons. We discuss it here for illustrative purposes. The second case is when the 
boxes are in a number eigenstate and the reservoir in a coherent state, and finally we 
consider when the boxes are in a coherent state as well as the reservoir. 



3.1. Superradiant Tunnelling Between Number Eigenstates 

The simplest case is when both the boxes and the reservoir are in eigenstates of the 
number operators , respectively. The unperturbed hamiltonian is: 

H = E^S^ + 2^FrS^ (18) 

The eigenstates of this hamiltonian are product states of the number states |A^r.) 
and \Nb). These states are equivalent to the usual quantum spin eigenstates \J,m) with 
J = 1/2 and m = N -1/2. 

This hamiltonian is perturbed by the tunnelling hamiltonian: 

Ht = -T{S^S; + SbS^) (19) 

and the current, i.e. the rate of change of charge on the island, is given by the 
commutator of the number operator with the hamiltonian: 
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dN^ _ dS^ _ 1 rr^ 

~dr~^~m'^''^^ 

i = T{sts; - s^st) (20) 

The expectation value of the current operator with the ground state is equal to 
zero, so we go to time dependent perturbation theory: 



m)) = u{t)m)) 
m)) 



i^ + ^J e-'''''/^HTe-''''''^dt' + • • -j |^/;(0)) (21) 
The expectation value of any operator at a given time is, to first order in Ht'- 

m)\im)) = {m\i\m) 

+ 23?{(^(0)|/4 / e'^*'/^i^Te-'^*'/Mt'|V^(0))} (22) 
in J 



The tunnelling hamiltonian and current contain only raising and lowering operators, 
and so their action on a number eigenstate produces a simple sum over two other number 
eigenstates, with a time dependent factor from the exponentials. 

The current is: 

-2T2 * 



mmm) = -r^ J cos (« - 2^^,) t'/h)dt' 

^ ' 



X ({N,\StS^\N,){Nr\S;St\Nr) - {N,\S^St\N,){Nr\StS;\Nr)) 

(23) 

The time dependence is determined by the energy difference between adjacent 
number eigenstates, i.e. by Ef^ and 2^Fr- If we evaluate the expectation values, we 
find terms like A^6(Z6 — Nij){lr — 2Nr) which is quadratic in If, and linear in Ij. (and vice 
versa, i.e. terms quadratic in Ir and linear in li,). In particular, to have the largest effect, 
we set Nf, — lf,/2 i.e. half the boxes are occupied. 



_2T2sin (EP-2^Fr)^A //, L \ 

^ w (!i;-2,l)/.) (i'l + - ^"'^ 

(24) 

We see that the current is proportional to l^, the square of the number of Cooper 
Pair boxes. This is in contrast to the situation where the initial state of the boxes is 
a product state of independent wavefunctions, in which any current can be at most 
linear in the number of boxes. The current is also proportional to (n^ — 1/2)/^, where 
Ur = Nr/lr, i-G. how far away the large superconductor is from half filling. 
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3.2. Superradiant Tunnelling Between a Number Eigenstate and a Coherent State 

The second case we consider is when the Cooper Pair boxes are in an eigenstate of the 
number operator, and the large superconductor is in a coherent, that is to say BCS-hke, 
state. We ensure the coherent state of the large superconductor by introducing a pairing 
parameter into the unperturbed hamiltonian: 

H = El'^S^ + 2iFrSr - - /^S; (25) 

The ground state of the large superconductor is a coherent state and = (5*7) 
can now be found self-consistently (see section ESI)- 

We make a basis transform and write the hamiltonian in terms of this new basis, 

h: 

H = El^S^ + 2EFr.S^^^ (26) 

This basis transform makes the unperturbed hamiltonian simpler, but at the 
expense of making the perturbation more complicated in the new basis: 



= — TSfj (— 5^ ,^ + fj. S + CrS fj) 

J^Fr 
2A* z 

J^Fr 

Again, the tunnelling current is zero to first order in T, as the boxes are in a number 
eigenstate. With this basis transform, we can easily calculate the second-order current. 
Following the same procedure as before, we get: 



immm) = -^^(o.,ni5.^^5.^jo.,.) j dfcos^f f/n 



X {{NhlS^S^m - {N,\S^S^\m)) 
2T2 



:{^r,h\S^f^S^fi\Or,h) 



{thy 
t 

{\e\'J dt' {cos{2EFr - Et) t'/h) {Nb\S+S^\Nk) 



t 

- \f\'J dt'cos((2E^, + Ef ) t'/n){m\S;S^\N,)} (28) 



The time dependence is again given by the level spacing, but for the large superconductor 



this is now given by Epr = y^Pr + I^P rather than {epr — A*)- 

If we assume the levels of the large superconductor are much more finely spaced 
than those of the boxes, i.e. Ep^ <^ E^^^, then: 
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mmm) (29) 

-2T^ sin (Eft A) |A^,p // \^ 

2T2 sin (EfH/h] , f„ , 

Remembering that Ir^pr/Ep^ = — N,., and setting A^^^ = hf^, we find this is 
proportional to: 

/.(l-2n,,)|(| + l) (30) 



This is very similar to the case described in ()24|) where both the boxes and the 
reservoir are in number eigenstates. In both cases, the effect is largest when the boxes 
are at half-filling, and the large superconductor is far from half-filling. The difference is 
that here it is the average number Nr that is relevant, not the number eigenvalue. 



3.3. Superradiant Tunnelling Between Coherent States 

The third case is where both the boxes and the large superconductors are in the coherent 
state. Pairing parameters are introduced for each: 



H — 2^FbS^ — AfeSf^ — A^Sft 

+ 2^FrS^ - ArS^ - A;S; (31) 

The ground state of this hamiltonian is for both superconductors to be in a coherent 
state, which again is determined self consistently. As before we make a change of basis, 
this time for both the boxes and the large superconductor, which makes the unperturbed 
hamiltonian: 

H = 2EFbS^f^ + 2EFrSr^f^ (32) 
The tunnelling hamiltonian becomes: 



Ht= - T{s+s; + Si;s+) (33) 

Al z Z 
J^Fb ^Fr 

Af) A* 

hipb ' ^Fr 

If we insert this into (|22p we find that we do have a non-zero term linear in T, i.e. 
the expectation value of / with 17/^(0)). To second order in T we have terms proportional 
to Ir and quadratic in and vice versa, and terms linear in both, i.e. 
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mmm) = {mimm + {m\i\m)iii. 

The first term is just tfie expectation value of / with the ground state: 
(^(0)|/|^(0)) = Ji^iM 2 sin (0, - 0,) 



ih 2Epfj lEpy. 
T lAJ |A 



ih Vh Vr 



2 sin — (j)r) (35) 



where 06, 0r are the phases of A;,, A^. This is the usual Josephson effect, as expected 
between two BCS-like states. To see how the tunnelling is enhanced by superradiance, 
we examine the addtional terms, i.e. the terms quadratic in and Ir- 
The term quadratic in is: 

mmm)iiw 

_ -T^ 2, Sm2EFrt/h\Ab\'^ ^pr 
LuLr 



{ihy ^ " 2Epr/h 2El^ Epr 
-T^ 2, cos2Eprt/n lAbplA, 



2Epjn 4Ei:4E|;'""'^'^^ - '^^^^ ('^^ 

The first term is proportional to |Abp/i?|.^ and ^pr/Epr, and so is largest when the 
boxes are half occupied and the large superconductor is far from half occupied. This 
term is also phase- independent and is like the effects seen in sections 13.11 and 13.21 

The second term is largest when both are half occupied, and also depends on the 
phase difference of the two pairing parameters with a frequency twice that of the usual 
Josephson effect. As such, it can be considered a higher harmonic in the Josephson 
current. The term quadratic in is identical apart from an anti-symmetric exchange of 
the labels, b, r. 

Finally, the term linear in both contains a phase - independent term that is largest 
when both the boxes and the reservoir are far from half-filling, and a term that is phase 
dependent with twice the Josephson frequency that is largest at half filling. 

{m\i\mkir (37) 

_ -T2 sm2{Epb + Epr)t/h 
~{in? 2{Epk + Epr)/h 

T2 cos2{Epb + Epr)t/h 

^ ' ' 2{Epb + Epr)/h 

In summary we conclude that tunnelling into or out of a Cooper Pair reservoir from 
a coherent ensemble of Cooper Pair Boxes can lead to a tunnelling current proportional 
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to the square of the number of 'boxes.' Interestingly, an observation of such scahng, 
in turn, could be taken as a demonstration of coherence and entanglement of the the 
'boxes' as such states are necessary to produce the phenomena. 



4. A 'Superconducting Jaynes-Cummings Model' and Revival 

4.I. The Superconducting Analogue of the Jaynes-Cummings Model 

Another well-known quantum optical phenomena is that of quantum revival (see e.g. 
[H]), which occurs when an atom is coupled to a single mode of the electromagnetic field. 
The initial state of the atom decays, but reappears at a later time. This apparent decay 
and subsequent revival occur because the environment of the atom is a discrete quantum 
field. Here we examine a Josephson system analogue. We consider a hamiltonian (eq. 

discussed in section \'A.2\ where the hamiltonian for the box has a number state as a 
ground state, and the reservoir ground state is a BCS-like state. These are the ground 
states we would expect for an uncoupled box and reservoir. The tunnelling term Ht 
couples these states. 



H = El'^S^ + 2iFrS^ - AS+ - A*S; 

= — T[Sj^ Sj. + Sf^ S^) (38) 

As in section 13. 2| we make a basis transform so that the reservoir is written as a 
quantum spin in another direction. 



zj rpch qZ I o 7-1 qZ 

Ht = — TSfj {— — S + f*S + CrS f^) 
2A* z 

~ "^^b (t; ^r,h + ^r^r,h + frS^A) (39) 

We can write a general state of the system in terms of the number eigenstates: 

m))= E a^,.^^MNr)\N,) (40) 

For a single two level system, = 1, but we can also do the calculations for a general 
number of symmetrised boxes. We can calculate the time evolution of the system using 
the Schrodinger equations: 



d a^^^N^it) 
Jt 



Y: a^^,r.St) (E^b'iNt - |) + 2MNr - |)) \Nr)\N, 
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OA 

Nr,Nf, 

OA* 



(41) 



More explicitly, we have a set of coupled differential equations for the coefficients 
aNr,Nti't)- 

^^da^^ = a^^Mt) [E^im -\)+ 2^MNr - |)) 

- a;v.,iv,+i(t) tJ{N, + l){h - iVfe)^(iV, - /,/2) 

- a^,_i,^,+i(t) Ty'(iV, + 1)(4 - N,)e;^jNrilr - iV, + 1) 

- a,v.+i,jv,+i(t) Tv'(iV, + 1)(/, - iV6)/,v'(iV, + 1)(/, - Nr) 

- a^.,iv,-i(t) T7iVb(4-iV,+ l)— ^(iV, - /,/2) 

- a^_i,jv,-i(t) T^iVfc(4 - iV, + l)f:^Nrilr -Nr + 1) 

- aN^+i,N,-i{t) T^N,{k -N, + l)er^{Nr + 1)(/, - iV,) 

(42) 

These equations are tri-diagonal and can be solved numerically. For now, to 
stimulate interest, we consider an approximation that allows a simple analytical model. 

The A terms in eq. EHlniean that the total particle number Nr+Nb is not conserved. 
We note that whilst the commutator of the total number operator [Nr + Nb, H] = 
AS^ — A*S~ is not zero, its expectation value in the coherent state is. With this in 
mind, we approximate and assume no fluctuations in the total number (although we of 
course retain fluctuations in N,. — Nb) and discard any terms in eq. 021 that change the 
total number. 

With these terms discarded, we have a set of coupled differential equations in which 
each coefficient a^^^N^it) only couples to a^^^i^N^^iit) and aNr+i,Ni-i(t). 

^^da^^ = a..,.,(t) [EtiNb - |) + EMNr - |)) 

- aN^+i,N,-l{t) T^Nbih -Nb + l)^{Nr + 1)(/. - Nr) 



- aNr-l,N,+l{t) T^{Nb + mb - Nb)^Nr{lr -Nr + l) 

(43) 

These equations can be restated in terms of an eigenvalue problem. Each level Nr 
has a set of eigenvectors for the boxes. If the boxes and the reservoir are initially in a 
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product state, the probability to be in a given value of Ni, (regardless of Nr) is: 

k 

Pn, = E I«^.(o)I' 

Where l^i^Nr) (Ei^Nr) are the eigenvectors (values) of the A^^-level system associated with 
the reservoir level \Nr), aAr^(O) are the initial amplitudes of \Nr), and I'lpNtiO)) is the 
initial state of the boxes. 

This may not be exactly solvable in general, but can be solved for few-level systems, 
lb- For example, if we have a single two level system and specify the initial state to be 
a product state of the box state \0b) with some state J2o.Nr{^)\Nr) of the reservoir, we 
have the probabilities that the box and reservoir are in a given state: 



-lEi 



(44) 



Pi,Nr-i{t) = |a;v,.(0)p-^sin2(n^,.t/^) 

Po,7v.(t) = |a^^(0)p(l-Pi,^,_i(t)) (45) 



With = TJNr{lr-Nr + l) and VLn^ = J {E^^ - Eprf / A + T'^ . This i 



IS 



the same result as in quantum optics, except that the photon creation/ annihilation 
operators give T'^ = Ty^I^ 



4-2. Quantum Revival of the Initial State 

In quantum optics, the phenomenon of revival is seen when the electromagnetic field is 
placed in a coherent state (i.e. la^l"^ = exp{—N)N^ /N\, with the sum over N running 
to infinity). We place the reservoir in the spin coherent state, i.e. 

Where a is determined by the average number |a| = Nr/{lr — Nr). 
The probability of the box being in the state |0b) at a given time is: 

P,{t)=J2\a^M\'PoMt) (47) 

For simplicity we consider the case when E^^ = Epr, i.e. fijv,, = T^v,.- This is shown 
in Figure El for the values Nr = 10, Ir = 50. The initial state dies away, to be revived 
at a later time. 

The revival occurs when the terms in the sum are in phase. We can make an 
estimate of this by requiring the terms close to Nr to be in phase: 



27r = 2T^{Nr + l){lr - Nr) trev " 2T^jNr{lr - Nr + 1) trev 

_27r ny^(l-n.)V^ 
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1 2 3 4 5 6 

time 



Figure 2. Spin Coherent State Revival for Nr = 10, Ir = 50 (solid line). The faint 
line indicates the usual quantum optical coherent state revival with A^^^ = 10, scaled 
along the t axis by and the vertical dash-dot line indicates the calculated revival 
time. 



Where, as before rir — Nr/lr- Note that the revival time remains finite in the 
Ir ^ oo hmit, as long as Nr/lr is finite. 

We can find an analytic form for the initial decay, and also demonstrate the 
equivalence of this 'spin revival' to the well known quantum optics revival, in the limit 
Ir — > oo, Ir/Nr — > oo. Note that this is an unrealistic limit for our case, where we would 
keep the filling factor Nr/lr constant, but it allows us to make a connection with the 
quantum optics case. 

The coefficients, aiv^(O) for the coherent and spin-coherent states are the Poisson 
and binomial distributions respectively, and both can be approximated by a Gaussian 
distribution. 



exp(- 



-Nr)Nl:'^/Nr\ - (27riV,)-i/2 



{2T:Nr) 



-1/2 



exp 



exp 



-{Nr - Nrf /2Nr 
■{Nr- Nrf/2Nr 



(49) 



For the coherent state we have taken the limit iV^ ^ 1. For the spin-coherent case 
we have also assumed Ir S> Nr- The gaussian factor suppresses any terms not around 
the average Nr, so we can expand Nr o::^ Nr + {^r — ^r) ■ 

In the limit, the spin coherent state gives: 



1 



Po{t) = o + 



1 



2 2(27rAr)V2 



dA^ exp 



{Nr N,)' 
2Nr 
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Figure 3. The large Nr,lr limit: Solid line indicates single box revival for Nr — 10, 
Ir — 100, and the dashed line indicates the asymptotic form for the decay. 



X COS 



2Nr 



(50) 



This is the same as for the coherent state apart from the factor of if the 
frequency. Doing the integral gives: 

Po(t) = 1 + 1 cos{2T (IrN^y/H) exp(-^^) (51) 
This is plotted in Figure El 

Thus the phenomenon of quantum revival has a direct analogy in our 'spin 
superconductor' model. In both cases, the revival is due to constructive interference 
between the terms in the sum over the reservoir (field) number states. 

4-3. Revival of Entanglement 

Revival is usually considered for single two-level systems, but the formula in (jl^ gives 
the state for a general system. In particular, the dynamics of a pair of two-level systems 
can be easily calculated, either numerically, or analytically, remembering that the singlet 
state is uncoupled due to radiation trapping. 
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Figure 4. Two Cooper Pair Box Spin Coherent State Revival for — 10, Ir — 50. 
The sohd line represents the probability that, when the system is initialised in the 
state |00), one of the Cooper pair boxes is in the state |0) at time t, after tracing over 
the other box. The dashed line indicates the revival of a single isolated Cooper Pair 
box connected to a reservoir. Note that the revivals occur at the same time, but out 
of phase. 

Starting in the state |00), we find a very similar effect to the single-island case. 
When one island is traced out, the revival of the other occurs at the same time as a 
single island, but with opposite phase (Figure IH). 

If we do not trace over one island, but consider the revival of the state |00), we 
find oscillations in the probability of returning to this state earlier than in the case of a 
single island. 

To investigate the revival of entanglement, we start the system in the state 
(|01) + |10))/v^- As before we note that the oscillations in the probability of being 
in the state (|01) + |10))/a/2 decay and then revive. It is interesting to note, however, 
that in this case the initial state is an entangled one. One measure of entanglement is 
the negativity i.e. the sum of the negative eigenvalues of the partially transposed 
density matrix. If we plot the negativity over time, we see that it dies away at first, as 
do the oscillations. As the oscillations revive, we also see an revival in the negativity of 
the two islands (Figure El). 

This revival of entanglement, an initially counterintuitive phenomenon, is of course 
due to the fact that we have considered the entanglement between the two islands only. 
When the entanglement 'disappears,' it is only the entanglement between the islands 
that has disappeared. If we were to consider the whole system, including the reservoir, 
the evolution would be unitary and the entanglement would remain constant. 

It is worth noting that this 'entanglement revival' is not a unique feature of our 
spin model, and would be expected whenever two-level systems are coupled to a field 
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Figure 5. Two Cooper Pair Box Spin Coherent State Revival for A',, = 10, Ir = 50. 
The dashed hne represents the probabihty that the Cooper Pair boxes have returned 
to their initial state (|01) + 1 10))/ ^/2 at time t. The dashed line indicates the negativity 
of the boxes as a function of time. 

with (a possibly infinite number of) discrete energy levels. 
5. Conclusions 

We have presented a simple model of a superconductor that allows an easy investigation 
of various phenomena associated with fluctuations around the BCS, mean field, 
behaviour. By taking the limit where all the single-electron energy levels are equal, 
we can represent the superconducting hamiltonian in terms of the operators of a large, 
finite quantum spin, which we then consider in the limit Ir ^ oo. Comparison of results 
using this model in the mean field with the limit of BCS results gives us confidence 
that whilst the model can be regarded as a caricature of a realistic description it retains 
some generic features of quantum fluctuations. Furthermore, we have demonstrated that 
the simplicity of the model allows the easy investigation of several surprising physical 
phenomena. 

As the first of these we considered an analogue of the quantum optical effect 
of superradiance, where the emitted intensity of a number of atoms coupled to an 
electromagnetic field is proportional to the square of the number of atoms. 

In the analogous ' Super- Josephson effect,' the current from a set of Cooper Pair 
boxes connected to a reservoir is proportional to the square of the number of boxes 
This effect is due to the entangled nature of the initial state. Any product state (or 
mixture of product states) would produce a current at most linear in The effect also 
relies on the reservoir being away from the Nr/lr = 1/2 state, i.e. away from half-filling. 

The analogy has been extended to the physically relevant case case where the 
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reservoir is in a coherent {i.e. BCS-like) state and the boxes are in a number state. As 
before we have no current to first order in T, but have a current proportional to to 
second order. Also as in the previous case, the effect requires the reservoir to be far 
from half-filling, but here it is the average number that must be far from lr/2. 

Finally we considered the effect when both the reservoir and the boxes are initially 
in a coherent state. In addition to the usual Josephson current, we observe the effect 
described above and also find a current that is dependent on the phase difference between 
the superconductors, with twice the phase dependence of the usual Josephson effect. 

Interestingly, the above discussion suggests that an observation of a Josephson 
current which scales quadratically with the number of small superconductors (qubits) 
can be taken as evidence for their entanglement. As experiments are now at the level 
of entangling two charge qubits [S], and these effects can be detected by making only 
current measurements, the above effects could be usefully looked for experimentally. 

We have also studied an analogue of quantum revival (again considering the case 
a box in a number state coupled to a reservoir in a BCS-like state) , in which the 
oscillations (of the type observed in [3111010]) in the probability of occupation of an 
island decay at first, only to revive at a later time. This effect can also be seen when 
we have two Cooper Pair boxes coupled to a reservoir. If we place these boxes in an 
entangled state and calculated the entanglement over time, we see that the entanglement 
between the boxes dies away at first, only to revive later. This remarkable effect 
demonstrates and measures the quantum fluctuations of the reservoir superconductor 
about the BCS mean field theory. In other words it probes the discreteness of the 
Cooper Pairs in the same way as the corresponding quantum optical effect probes the 
discreteness of photons of the radiation field. 
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Appendix A. Limit of BCS Self Consistency Equations 

As confirmation of the validity of our spin model, we wish to compare it with the 
standard BCS theory. In particular, we wish to check that solving the gap equations 
from the spin model gives the same result (fT^ IT!^ as if we solve the usual BCS gap 
equations first, and then take the limit. We wish to take the limit where the single- 
electron energy levels, are equal. To do this we consider a density of states which is 
a top-hat function centred around ep with width uj and height Ir/uj. Taking the limit 
uj —>■ ensures that the interaction occurs in a narrow region around ep whilst ensuring 
the integral is over a constant number of levels, Ir. We have two integral equations that 
must be solved self-consistently. The first is an equation for the pairing parameter. A: 
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and the second is an equation for the average number of Cooper Pairs, N: 

iV = i / de 1 '-^ ,] (A.2) 

As we are at zero temperature, both these integrals can be done exactly: 

(A.3) 



1 


I 








In 


V ~ 












N = 


-( 











(eF + ^/2) + v'(eF + ^/2)2+|A|^ 



{^F + ujl2f + |AP + ^(e^ - ujl2y + |A|2j (A.4) 

Again, is shorthand for ep — fJ'- These can be rearranged to give two equations 
for |Ap ffX5irO|l . or two equations for ^j. ffX8irO|) . The equations for |Ap are: 



Equating the two and rearranging gives us an expression for ^p'^: 



where x = evi and n = —. If we rearrange eqns. ()A.3l IA.4|) for we get: 



(l-2„)^ + |A|^A__^ (A.9) 
Which yield: 

From eqns. ()A.7l IA.10|) . it is clear that the correct limit to take is u/V 0. Taking 
this limit of eqns. (f^rOn|l . we get: 

T/2/2 

=^(l-2n)^ (A.ll) 

|A|2 = (A.12) 
That is, we have regained the forms for A and ^ given by the spin model (jl2t fT ^ . 
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